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We give a detailed microscopic analysis of why holes are different from electrons in condensed 
matter. Starting from a single atom with zero, one and two electrons, we show that the spectral 
functions for electrons and for holes are qualitatively different because of electron-electron interac- 
tions. The quantitative importance of this difference increases as the charge of the nucleus decreases. 
Extrapolating our atomic analysis to the solid, we discuss the expected differences in the single par- 
ticle spectral function and in frequency dependent transport properties for solids with nearly empty 
and nearly full electronic energy bands. We discuss the expected dependence of these quantities on 
doping, and the physics of superconductivity that results. We also discuss how these features of the 
atomic physics can be modeled by a variety of model Hamiltonians. 



I. INTRODUCTION 

The understanding of electronic correlation in solids 
and its consequences for charge transport and for collec- 
tive phenomena are fundamental problems in condensed 
matter physics. At the simplest level, electronic correla- 
tion manifests itself in the two-electron atom [Q . In this 
paper we return to the old problem of the two-electron 
atom Q and show that a fundamental physical principle 
emerges from it, electron — hole asymmetry. We then 
show that the understanding of electron-hole asymmetry 
at the atomic level provides insight that implies that this 
asymmetry has fundamental consequences for the physics 
of charge transport and collective phenomena of the solid 
state. 

Remarkably, none of the many-body model Hamilto- 
nians that are most widely used to study the effect of 
electronic correlation in solids such as the Hubbard 
model, extended Hubbard model, degenerate Hubbard 
model, the Anderson impurity and lattice models, the 
Kondo model, Falicov-Kimball model, Holstein model, 
Su-Schrieffer-Heeger model and t — J model, contain this 
very basic and fundamental aspect of electronic corre- 
lation that follows from the atomic analysis. Hence we 
argue that these models are fundamentally flawed. We 
propose a variety of other model Hamiltonians that con- 
tain this physics, which may be generically called 'dy- 
namic Hubbard models' |5|-||. These models lead to a 
new understanding of the physics of charge transport in 
solids, and to a universal mechanism of superconductiv- 

ity §0. 

The observation that electrons and holes in atomic 
shells are in some sense equivalent was first made by 
Heisenberg |]l2| . It is easy to see that in the absence of 
electron-electron interactions, an atom with i electrons in 
an outer partially filled shell defined by quantum num- 



bers n, I, has the same multiplet structure as an atom 
with p — i electrons in this shell, with p = 41 + 2 the total 
number of electrons that can fill this shell. In particular, 
this is trivially true for an s-shcll (1=0), where the state 
with one electron in the ns° shell is identical to the state 
with one hole in the ns 2 shell. Furthermore, all the ma- 
trix elements of the Coulomb interaction operator, both 
diagonal and off-diagonal, between these non-interacting 
electron atomic states, are the same for a shell with i elec- 
trons and one with i holes (i.e. p — i electrons) |Q. This 
of course extends also to cases where more than one shell 
is incomplete. As a consequence, the atomic multiplet 
structure of a shell with i electrons is identical to the one 
with i holes within first order perturbation theory in the 
Coulomb interaction between electrons. Thus, to the ex- 
tent that first order perturbation theory holds for atoms, 
electron- hole symmetry exists for the atomic states, and 
similar arguments can be used to argue for electron-hole 
symmetry in energy bands in the solid state jj3. This 
implies that the states of atomic shells that are almost 
full, and of electronic energy bands that are almost full, 
can be described in terms of the states of few positively 
charged holes, instead of the states of many negatively 
charged electrons. The physical reality of the hole con- 
cept for transport in solids was first demonstrated by 
Hall's measurement |ll| and Peierls' explanation [|l6| of 
positive Hall coefficient in solids with nearly filled elec- 
tronic energy bands. 

However, in this paper we show that electron-hole 
symmetry in fact does not exist in atoms, and conse- 
quently it will not exist in electronic energy bands in 
solids. Instead, we argue that electrons and holes are 
fundamentally different objects. In some sense, in the 
light of the above discussion, this can be understood as 
a failure of first order perturbation theory. It is of course 
not surprising that perturbation theory in the electron- 
electron Coulomb interaction should fail. The electronic 
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stucture of non-interacting electrons in atoms is deter- 
mined by the Coulomb interaction between electrons and 
ions, which of course has the same coupling constant as 
the electron-electron Coulomb interaction. Thus, it is 
never true that the spacing between non-interacting elec- 
tron atomic energy levels is larger than matrix elements 
of the electron-electron Coulomb interaction, the regime 
where perturbation theory may be expected to be valid. 
Despite of this, it is fortunate that the qualitative struc- 
ture of atomic multiplets and their quantum numbers are 
the same for electron and holes |Q. 

The physics of electron-hole asymmetry that results 
from electron-electron interactions occurs even for s- 
shells, and even for the atomic Is shell, that is most 
isolated in energy. That is the case that we will analyze 
in this paper. It is natural to expect that the physics 
that we find here should have an even larger effect for 
atomic states with higher quantum numbers, where shell 
occupation becomes larger and energy level separation 
becomes smaller. 

Since for an s-shell the one-electron state and the one- 
hole state are the very same state, one may ask how it is 
possible to understand electron-hole asymmetry based on 
it. The point is, in considering the process of transport 
of charge in solids with electrons or holes, we need to 
consider not only the atom in a given charge state, but 
also processes where electrons and holes are created and 
destroyed as conduction occurs. It is these creation and 
destruction processes that contain the physics of electron- 
hole asymmetry even at the level of a single atom and 
even for s-shells. 

This paper is organized as follows. In the next sec- 
tion we discuss the electronic state of two electrons in 
the Is shell of ions of nuclear charge Z . We show that 
the state of the two-electron atom is strongly influenced 
by electronic correlation effects, the more so the smaller 
Z is. In section III we discuss the calculation of the 
atomic spectral function for electrons and holes, and the 
qualitative reasons for why they are fundamentally dif- 
ferent. Section IV discusses results for the atomic spec- 
tral function within the various approximations for the 
state of the two-electron atom reviewed in Sect. II. In 
Sect. V we discuss the consequences of these results for 
the properties of electrons and holes in the solid state, 
in particular for the single electron spectral function and 
for frequency-dependent transport. In section VI we dis- 
cuss the consequences of this physics for the understand- 
ing of superconductivity in solids. Sect. VII discusses 
several model Hamiltonians that contain the physics of 
electron-hole asymmetry found at the atomic level. We 
conclude in Sect. VIII with a summary and discussion 
of our results, a review of the empirical and experimen- 
tal evidence in support of this physics as the underlying 
universal mechanism of superconductivity in solids, and 
a survey of some of the many open questions and oppor- 
tunities for further research in this area. 



II. THE TWO-ELECTRON ATOM 

The wavefunction for an electron in the lowest energy 
state (Is) of a hydrogen-like ion of nuclear charge Z is 



Mr) = ( — ) 1/2 e- Zr 
ir 



(1) 



with r measured in units of the Bohr radius ao. The 
ground state wave function of the two-electron ion is not, 
of course, 



(2) 



because of electron-electron interactions. Consider the 
following approximations to the ground state of the two- 
electron ion: 



A. The Hartree wavefunction 

The simplest approximate wavefunction that takes into 
account the effect of electron-electron interaction is of the 
Hartree form 



where 



#ff(ri,r 2 ) = ¥ J z(n)^(r 2 ) 
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(3a) 



(3b) 



results from minimization of the total energy for a vari- 
ational wave function of the form Eq. (3a). Eq. (3) 
describes an expanded orbital for each electron, due to 
the partial shielding of the positive nuclear charge by the 
other electron. In Eq. (3), both electrons reside in the 
same expanded orbital. A slightly better approximation 
of the Hartree form is obtained by allowing for the most 
general single particle wave function rather than the ex- 
ponential form Eq. (1). For example, the error in the 
ground state energy of the He atom is 1.95% with Eq. 
(3) and 1.45% with the optimal single particle wavefunc- 
tion in the Hartree wavefunction ffl]. 



B. The Eckart wavefunction 

A better approximation is obtained by allowing for ra- 
dial correlations, through the wave function jnj 



2(1 + ^)1/2 



S12 = (¥>Zi,Vz 2 ) 



(4b) 



where the exponents Z\ and Z2 are again obtained by 
minimization of the total energy, for a wavefunction of 
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the form Eq. (4). The numerical results for Z\ and Z 2 
obey the approximate relations 



Z x = l.UZ- 0.105 



Z 2 = 0.905Z- 0.622. 



(5a) 



(5b) 



Hence, one of the electrons resides in an orbit of approxi- 
mately the same radius as that of the single electron ion, 
and the second one resides in a substantially enlarged or- 
bit. The minimization procedure resulting in the values 
of Z\ and Z 2 becomes unstable for Z < 0.93. For He, 
the error in the ground state energy is now reduced to 
0.98%. 



C. The Hylleraas wavefunction 

Much more accurate wavefunctions for the two- 
electron system are obtained by introducing dependence 
of the wavefunction on r\ 2 = \r\ — 7^*2 1 in addition to r\ 
and r 2 , which allows for angular in addition to radial cor- 
relations. We consider here the simplest wavefunction of 
this form, Hylleraas' "dritte Naherung" (third approxi- 
mation) Q): 

^Hy(ri,r 2 ) = N(p(2Zks, 2Zkt, 2Zku) (6a) 



(f(s,t,u) = e~ s/2 [l + c 1 u + c 2 t 2 } (6b) 



s = ri + r 2 (6c) 
t = r 2 — r\ (6d) 

u = r 12 (6e) 

The parameters c\ , c 2 and k are determined by minimiza- 
tion of the energy, as described in the appendix, and N 
is a normalization constant. Both this as well as the pre- 
vious wave functions yield upper bounds to the ground 
state energy. C\ and c 2 are found to be positive, indicat- 
ing that for given s = r\ + r 2 the wavefunction is larger 
for large angle between f*i and r 2 (angular correlations) 
and for large \r\ — r 2 \ (radial correlations). The error 
in the He energy with the wavefunction Eq. (6) is re- 
duced to 0.05%. Table 1 summarizes the values of the 
energy and orbital exponents for the three wavefunctions 
discussed for the He (Z = 2) and the H~ (Z = 1) ions 
compared to the experimental values. 

Figure 1 shows the orbital exponents versus ionic 
charge in the different approximations. The Hylleraas 
wave function has an orbital exponent similar to the 
Hartree wave function, while the exponents of the Eckart 
wave functions are very different. This suggests that even 
though the Eckart wave function gives better results than 



the Hartree one for the energy, the Hartree method may 
give a better representation of the wavefunction itself. 
Figure 2 shows the parameters c% and c 2 of the Hylleraas 
wavefunction, which describe angular and radial corre- 
lations respectively, versus ionic charge. Note that both 
parameters increase rapidly as the ionic charge decreases, 
and that angular correlations are much more important 
than radial correlations for large Z. 

The different wavefunctions discussed above describe 
attempts of the electronic wavefunction to reduce the 
Coulomb repulsion between electrons, without paying un- 
duly in electron-ion energy. The effective Coulomb repul- 
sion between electrons is 



U eff {Z)=E(2)+E(0)-2E(l) 



(7) 



with E(n) the ground state energy for the ion with n elec- 
trons (E(0) = 0, E{1) = -Z 2 , with E in Ry = 13.6eV). 
For the wavefunction Eq. (2) with the electrons occupy- 
ing the same orbitals as in the singly occupied atom 



Ueff(Z) = U=-Z 



(8a) 



U 



J d VdV|^(r)| 2 ^^|^(r')| 2 (8b) 



and for the Hartree wave function Eq. (3) 



25 
128' 



(9) 



smaller than Eq. (8a) because the wavefunctions are 
more expanded (note however that this is different from 
the actual Coulomb integral Eq. (8b) computed with 
the Hartree wavefunctions, due to the cost in electron- 
ion energy payed by the expanded wave functions). For 
the Eckart and Hylleraas wavefunctions the effective U 
becomes progressively smaller, corresponding to the de- 
crease of the two-electron energy E(2) with increasingly 
better wavefunctions. Figure 3 shows U e ff versus ionic 
charge in the different approximations. 

Figure 4 depicts qualitatively the electrons in the atom 
in the different approximations. We emphasize that the 
two-electron wavefunction in any of the approximations 
discussed is very different from the one corresponding 
to non-interacting electrons, Eq. (2). For the purposes 
of this paper, any of these wavefunctions, including the 
simple Hartree one, describes 'electronic correlation', in 
the sense that the two-electron wavefunction is different 
from the one given in Eq. (2), and is appropriate to 
illustrate the physics of electron-hole asymmetry. Note 
also that the difference between these wave functions and 
the uncorrelated wavefunction Eq. (2) increases as the 
ionic charge Z decreases and the relative importance of 
electron-electron versus electron-ion interaction strength 
increases. 

It is also worth pointing out that the electron-electron 
repulsion U e ff decreases as the ionic charge Z decreases 
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for all these wavefunctions, including the uncorrelated 
one. Hence, everything else being equal, pairing of elec- 
trons (or holes) should be favored when the effective ionic 
charge is small, i.e. with negatively charged ions. This 
simple fact, which provides some rationale for high tem- 
perature superconductivity being favored in systems with 
negative ions (0 = in cuprates, B~ in MgB-i) has surpris- 
ingly not been pointed out before to our knowledge. 



III. ATOMIC SPECTRAL FUNCTIONS 

Properties of many-electron systems can be studied by 
consideration of spectral functions. The zero tempera- 
ture one-electron spectral function for a many-body sys- 
tem is conventionally defined as 



A(u) 



OO 

D 

n=l 

I < rale, 



< n|cL|l > \ 2 5{" - « +1 < +1 + u N )) 
\1>\ 2 6(uj + (E»- 1 -EF- 1 -h N - 1 ))}. 



where 1 1 > is the ground state of the ./V-electron system, 



Ei its ground state energy, and /ijv = E^ +1 — E^ , 



Hn-i = E^ — E^ 1 . For a metal, /in — Mw-i = Mi 



and one can redefine the frequency u — > u> 
Eq. (10a) becomes 



u so that 



A(lu) = zS(uj - (e - /it)) + 4 mc (w) 



(11) 



where the first term describes the quasiparticle, with 
quasiparticle weight z, < z < 1, and the second 
term describes a continuum of incoherent excitations at 
higher energies. For a small system like an atom how- 
ever, the spectral function will consist of only discrete 
5-functions. Nevertheless, we can identify the lowest- 
energy 5-function as corresponding to the quasiparticle, 
and its coefficient as the quasiparticle weight. 



The fundamental asymmetry between electrons and 
holes follows immediately from consideration of these 
spectral functions. Consider the spectral function for cre- 
ating an electron in the empty atom, N — 0. From Eq. 
(10b), it is simply 



A ei {uj) = 5{uj) 



(12) 



so that it is totally coherent, with quasiparticle weight 
z = 1. 

Consider next the spectral function for creating a hole 
in the full Is shell: 



A h {uj) 



<n\c aA \l > \*8(w + E*-El) (13a) 



The excitation energies in Eq. (13a) are the ones of the 
singly-occupied atom, 



(10a) 

and the matrix elements arc 



A(u) = £[| < n l4r|l > ~ (E^ +1 < +1 )) (10b) 

ri=l 

+ | < n\c aa \l > \ 2 5(lu + «" x - E?- 1 ))]. 

and the zeros of the 5-function arguments occur at the 
excitation energies of the system. To have a smooth tran- 
sition between the metal and atomic spectral functions 
we will here define the atomic spectral functions as Eq. 
(10b) rather than Eq. (10a). The first term describes 
the response of the system when an electron of spin a 
is created in the single-particle state a, into the ground 
state with N electrons. E^ +1 is the n-th excited state 
energy of the many-electron system with N + 1 electrons. 
Similarly, the second term describes the response of the 
system when an electron of spin a at the single-particle 
state a is destroyed from the ground state of the N- 
electron system, and E^ 1 are the excited state energies 
of the resulting (N ~ l)-electron system. Generically, in 
many-body systems the one-electron spectral function is 
of the form M 



7 2 

E n — 2 



(13b) 



< H<Vt|l >= s n = J d 3 rid 3 r2^(fi,f2)(p a (ri)(p n (r 2 ) 

(13c) 

with ^ the two-electron ground state wave function, ip n 
the wavefunction for the n-th excited state of the singly 
occupied atom, and tp a the single particle wave function 
of state a. The important qualitative point is, because 
^ is not given by the product form of single electron 
wavefunctions Eq. (2), the spectral function Eq. (13a) 
necessarily consists of a sum of (5-functions rather than 
a single one as Eq. (12): 



A h (u) = zS(lu) + V S 2 J(u + E\- E\) 



n^l 



si 



(14a) 



(14b) 



Hence, the spectral function for hole creation Eq. (14a) 
is qualitatively different from the one for electron cre- 
ation Eq. (12). 

The foregoing remarks are true for any single-particle 
wave function (p a . Which is the appropriate (p a to use 
to define the hole spectral function? We argue that it is 
reasonable to use the single particle wave function that 
maximizes the quasiparticle weight of the hole spectral 
function, Sf (such a hole wavefunction will have minimal 
kinetic energy in the solid) . For the Hartree wavefunction 
Eq. (3a) it is simply given by 



and 



V<x{r) = ifz(r) 



Si= d 3 rip 1 (r)ip 2 {r) = 



{zzyi 2 
(W 



(15a) 



(15b) 
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with Z = Z — 5/16. For the Eckart and Hylleraas wave 
functions, we use a single particle wave function of the 
form Eq. (15a), with Z determined by maximization of 
the function 



(16) 



Si(Z)= / d- 5 r 1 d 6 r2^{ri,f 2 )ipz(ri)ipi(r 2 ) 



Note that the hole spectral function as defined by Eq. 
(13) satisfies 

J du;A h (w) =J2 S n = J d ^l[ J d?r 2 ^(r u r 2 ) Va (r 2 )f 

(17) 



In the Hartree approximation 



duAh(w) 



(18) 



while with the more accurate wave functions, the fre- 
quency integral of the hole spectral function is less than 
unity. This is because there is also a non-zero probabil- 
ity of creating an electron in the single particle state (p a 
into the two-electron atom. However that probability is 
found to be small. 

Finally, the atomic spectral function for the singly oc- 
cupied atom is of interest. We consider separately the 
positive and negative frequency parts. The spectral func- 
tion for destruction of an electron in the singly-occupied 
atom is: 

AM =J2\<n\c a ,\l> \ 2 5(u; + E° n E?) (19) 



which is of course also given by 



(20) 



and the spectral function for destruction of a hole in the 
singly occupied atom, or equivalently creation of an elec- 
tron in the singly occupied atom, is : 



A' h (u>) = '£(S' n ) 2 5(u>-(E 2 n -E 2 1 )) 



(21) 



where E 2 are the excitation energies of the doubly occu- 
pied atom, and 



5,, 



d 3 rid 3 r 2 ip a (ri)^(r 2 )* n (n, r 2 ) 



(22) 



with ^f n the wavefunctions for the excited states of the 
two-electron atom. Once again, the spectral function A' h 
is qualitatively different from A' ell as it contains an inco- 
herent part. The spectral function for the singly occupied 
atom as defined by Eq. (10b) will be either Eq. (20) or 
Eq. (21) depending on the value of a in Eq. (10b) and 
of the spin of the electron in the singly occupied atom . 

To estimate S' n we need the excited state wavefunc- 
tions for the doubly occupied ion. The excited states are 



approximately described by one electron in the ground 
state and the other electron in an excited state of the ion 
with orbital exponent Z, <p nZ , with Z given by Slater's 
rules lit: 



Z- 0.85, n= 2 



(23a) 



Z = Z-l,n = 3,4, ... (23b) 
The matrix elements are then given by (for n > 1): 

S'n = J d s rip a (r)>p n z(r) 
and the excitation energies by 

Z 2 



(24) 



Ei = -Z l 



(25) 



IV. 



RESULTS FOR THE ATOMIC SPECTRAL 
FUNCTIONS 



As discussed in the previous section, we calculate the 
hole spectral function for that single particle wave func- 
tion <p a that maximizes the quasiparticle weight, which 
we will call the 'hole wavefunction'. Figure 5 shows the 
orbital exponent of that state, a, versus Z. For the 
Hartree approximation a — Z, for the Eckart and Hyller- 
aas approximations a is somewhat smaller, i.e. the hole 
wavefunction is somewhat more extended. In all approx- 
imations the hole wavefunction becomes more diffuse as 
the ionic charge Z decreases, as expected. Figure 6 shows 
the quasiparticle weight for the hole, from Eqs. (14b) and 
(16). The Hartree wavefunction somewhat overestimates 
and the Eckart wavefunction somewhat underestimates 
the quasiparticle weight given by the Hylleraas wavefunc- 
tion, which presumably is the most accurate one. As the 
ionic charge decreases the quasiparticle weight for the 
hole decreases in all the approximations. 

To obtain the full spectral function for hole creation we 
compute the matrix elements S n , Eq. (13c). Because ip a 
is spherically symmetric, only the excited atomic s-states 
give non-zero results for the integral. They are given by 



<Pn(r) 



1=0 



a, 



{n-l)\2 l 



(n-Z -!)!(/ + 1)!J! 



(26a) 



(26b) 



The integrals for all the approximate wavefunctions con- 
sidered can be done analytically, as described in Ap- 
pendix B. 

Figure 7 shows the spectral function for hole creation in 
the full Is shell, Eq. (13a), for various values of Z. We 



5 



show results for both the Hartree and Hylleraas wave- 
functions; the results for the Eckart wave function are 
similar. For Z — 2, Ah(uS) is close to a 6— function at 
lu = 0, of weight unity; the incoherent part occurs at very 
high energies and has very small weight. As Z decreases, 
the weight of the w — peak (quasiparticle) gradually 
decreases, and the weight of the higher enery peaks (in- 
coherent part) gradually increases; in addition, the inco- 
herent peaks shift to lower frequency. In Fig. 8 we show 
results with the Hartree wavefunction for an even smaller 
Z, Z = 0.4 (the procedure to obtain the Hylleraas wave- 
function does not converge for such a small Z), showing 
that here the incoherent contribution is bigger than the 
quasiparticle part. 

Similarly we can calculate the spectral function for hole 
destruction (electron creation) in the ion with one hole 
(i.e. with one electron), from Eqs. (23)-(25). The ex- 
cited states wavefunctions are given by Eq. (26) with Z 
replaced by Z given in Eq. (23). The matrix elements 
turn out to be similar to the case of hole creation, but 
the excitation energies here are much lower. Hence the 
incoherent part of the spectral function is shifted to much 
lower energies for these processes. Figure 9 shows the cal- 
culated spectral function for hole destruction for Z = 2 
and Z = 1. 

Note that it is meaningless to calculate the spectral 
function for hole destruction for the isolated atom for 
Z < 1. In that case, all the excited states of the two- 
electron ion are unbound. However the situation is dif- 
ferent in the solid state. There, the negative ion (e.g. 
= or B~) will be surrounded by positive cations, and 
excited states can be formed where the excited electron 
is still confined in the neighborhood of the anion and the 
surrounding cations. 



V. CONSEQUENCES FOR THE SOLID STATE 

We have seen in the last section that the spectral func- 
tions for electrons and holes in the atom will be either 
single (5-functions or sums of 5-functions depending on 
the charge state of the atom. In this section we discuss 
the fundamental consequences of this for the solid state. 

When we bring together the atoms to form a solid, 
orbitals overlap and bands are formed. Consider for sim- 
plicity the lowest band, formed by overlaping Is orbitals, 
well separated from other bands if the interatomic dis- 
tance is large. The kinetic energy operator for an electron 
in such a band, in second quantized formalism, is 



and c\ a creates an electron in Wannier orbital i, a linear 
combination of nearby atomic orbitals with largest am- 
plitude at site i. In the simplest tight binding approxi- 
mation the hopping matrix element for a single electron 
can be obtained from 



-(<Pi, hipj) 



Sij(<p,h<fi) 



q2 



with 



h = -V 2 - 



2Z 



2Z 



Rj\ 



(28a) 



(28b) 



(29) 



the single electron Hamiltonian (in atomic units). Eq. 
(28) reproduces quite accurately the spacing between 
bonding and antibonding states of the single electron 
diatomic molecular ion in a wide range of interatomic 
distances and ionic charges []l9| . However we emphasize 
that Eq. (27) is generally valid for any Bloch band: away 
from the tight binding limit the hopping will involve 
many neighboring atoms and differ from Eq. (28), and 
the Wannier orbital will be a linear combination of many 
atomic orbitals. We will focus on the tight binding limit 
here from simplicity, but expect the qualitative physics 
to survive beyond that limit. 



A. Spectral function 

1. General considerations 

Consider the spectral function for electrons in this 
band, A(k,uj). For the empty band, we use Eq. (10) 
with c Qcr replaced by 



Cfecr = 



—y 



(30) 



When c\ a acts on the ground state (the empty band), it 
yields a single electron Bloch function, an eigenstate of 
the Hamiltonian Eq. (27a). Hence the spectral function 
is simply given by 



A(k,uj) = S(uj - (e k - fx)) 



(27a) 



a- 



J_y^ ik(R i -R j ) t . 



(31a) 



(31b) 



where tij is the Fourier transform of the Bloch band en- 
ergy 



(27b) 



with fi the chemical potential (=the bottom of the band 
for the empty band). 

For the full band instead, the situation is much more 
complicated. We wish to compute Eq. (10) with the 
operator 



G 



c ko 



(32) 



where the atomic destruction operator cf a destroys an 
electron in the single electron atomic state a with wave- 
function cp a as discussed in the previous section. The 
ground state |0 > is given by 



|o>= i n>i i u> 2 ...| u>< ...| u> 



iV 



(33) 



where | TT>j is the correlated two-electron ground state 
of the i-th atom, described e.g. by the Hylleraas wave- 
function. When cf a operates on this state we obtain 



U>;=Si| i> l +J2 s n\ i n >i 



(34) 



n>l 



where the matrix elements S n are given by Eq. (13c) 
and (B7), and | [ n >i denotes the n-th excited state of 
the single electron ion with energy given by Eq. (13b), 
and | |>= | J. 1 >. Hence, 

c -i° >= 4E e ^i i>* + 7^E^ E s «i t n >* 

v i v i n>l 

(35) 

which is, not surprisingly, not an eigenstate of the Hamil- 
tonian Eq. (27a). (In Eq. (35), the state vectors for 
all sites j ^ i, \ ~\i>j, are omitted for clarity.) Thus 
the spectral function for a single hole in the full band is 
not a single 6— function and hence is qualitatively differ- 
ent from the spectral function of a single electron in the 
empty band, similarly to the case for the single atom. 

The first term in Eq. (35) denotes a Bloch wave for 
a single hole in the full band. This is also not an exact 
eigenstate of the Hamiltonian Eq. (27a), because every 
time the hole hops to a neighboring site there is a finite 
probability amplitude that the final state has the atoms 
in excited states, as discussed in the previous section. 
Nevertheless, within the approximation where those pro- 
cesses are excluded, the energy of the itinerant hole will 
be given by 



h o2 ~ 

-k — b l e k = z Ck — £fe- 



(36) 



Separating this contribution from the rest we have for 
the single hole spectral function 



Ah(k,w) = z5(u) - (ejfc - //)) + A'(k,uj) 



(37) 



where A'(k,ui) contains all contributions to the spectral 
function with atomic excited states resulting from the 
second term in Eq. (35). Hence it will involve the high 
energy part of the atomic spectral function discussed in 
the previous sections. In a solid, we expect that this con- 
tribution will not have a strong k-dependence and that 
the atomic ^-functions will broaden to give rise to the 
'incoherent' part of the spectral function, A'(k,u>). 



The derivation of Eq. (37) is not rigurous, and hence 
the expressions for z and are not exact. Neverthe- 
less, the general form Eq. (37), possibly with some k- 
dependence to the quasiparticle weight z, is expected to 
be correct for a many-body system [Q. In fact, Sham 
has shown rigurously (2(]] that the spectral function for 
a sigle hole in a full band has a ^-function quasiparticle 
contribution that is separated from the continuum of in- 
coherent excitations. We believe that the k-independent 
z and the form of obtained are reasonable approxima- 
tions at least in the tight binding limit. In particular, this 
implies that the quasiparticle weight z and the effective 
mass renormahzation for a hole 

- ^ - (38) 
m z 

will have the qualitative dependence on the ionic charge 
Z discussed in the previous section. 



2. Quasiparticle operators 

In the foregoing discussion we considered the spectral 
function for a single electron in the empty band and for 
a single hole in the full band. To generalize this analysis 
to other band fillings we define quasiparticle operators. 
The atomic electron operator can be written in terms of 
atomic states as 



Ci T = |o ><T | +SiU><u | + 



E 

(n,n')#(0,0) 



Snn'\ r><tr 



(39) 



where | |™> and | ||™ > are excited states of the singly 
and doubly occupied ion respectively, and 



S nn > = Id rd r ^n'{ri 1 r 2 )f a {r 1 )tp n {r2) 



(40) 



with S n i — S n , Sin' — S' n , as discussed in the previous 
section. We write Eq. (39) as 



('in 



[1 + (Si - tyfii-vjciv + c' ia 



(41) 



where the 'quasiparticle operator' Ci a acts as if electronic 
correlations didn't exist: 



4lo>= k> 



(42a) 



c l a\-<J>=\U> (42b) 

and obeys usual anticommutation relations. n.; CT = 
c ia .Cia-. The operator c' ia is the 'incoherent part' of the 
bare fermion operator Ci a , given by the last term in Eq. 
(39). From Eq. (41), we obtain that the quasiparticle 
weight for a band with band filling of n electrons per 
atom is approximately given by 



z{n) = [l + {S 1 -l)-f 



(43) 



decreasing monotonically from 1 to Sf as n increases from 
to 2. 
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3. Quasiparticle Hamiltonian 

The effective low energy Hamiltonian for quasiparti- 
cles (kinetic energy part only) is obtained by replacing 
the bare fermion operators in Eq. (27a) in terms of 
their form Eq. (41) and neglecting the 'incoherent part' 
of the operators that would leave the atoms in excited 
states. Hence, it describes only ground-state to ground- 
state transitions (diagonal transitions in the language of 
small polaron theory) 



H, 



qp 



(44a) 



excitation energies of the singly occupied ion and the 
doubly occupied ion, hence t\ > e 2 , and we assume the 
broadening I\ is proportional to e,. The relative weights 
of Ah and A' h in A' will be determined by the probabil- 
ities of doubly and singly ocuppied sites respectively in 
the solid for band filling n. The probabilities for empty 
singly and doubly occupied sites p e , p s , pd can be written 
as 



p e (n) = (1 - -f 



I n \2 

{ 2 )a 



p s (n) = n(l 



2 ) + 2a( 2 } 



(48a) 



(48b) 



tij[l - (1 - S , i)(ni,_ <T + 

+ (1 - Siffli-aflj-a]. 



(44b) 



In a mean field approximation, the quasiparticle energy 
is simply given by taking the expectation value of Eq. 
(44b), yielding 



e k = z{n)e k . 



(45) 



To simplify the notation, we omit to denote explicitely 
the n-dependence of e k on the left side of Eq. (45). Eq. 
(45) implies that the relation Eq. (38) between effective 
mass and quasiparticle weight is valid for any band fill- 
ing, and that the quasiparticle mass increases and the 
effective bandwidth decreases with increasing band fill- 
ing. 



4- Phenomenological model for the spectral function 

From the considerations in the previous sections we 
conclude that a reasonable model for the spectral func- 
tion in the solid that contains the physics discussed is of 
the form 



A(k, u) = z{n)5(uo - (e k - fi)) + A'(uj) 



(46) 



where A'(lj) is the incoherent part describing the atomic 
excitations, and we have ignored its k-dependence. The 
quasiparticle weight z(n) and energy e k are given by Eqs. 
(43) and (45) respectively. The atomic spectral functions 
discussed in Sect. Ill will be broadened in the solid state, 
and we assume the Gaussian forms 



A h (co) = 



4» 



1 



27iTi 



771— 



27rr 2 



(47a) 



(47b) 



corresponding to hole creation (negative u>) and hole de- 
struction (positive uj) respectively e x and e 2 are mean 



p d (n) = (-f{l-a) 



for n < 1, and 



p e {n) =pd(2- n) 



p s {n) =p 3 (2-n) 



Pd(n) =p e (2-n) 



(48c) 



(49a) 



(49b) 



(49c) 



for n > 1. Here, a is a parameter that describes the 
suppression of double occupancy due to the on-site re- 
pulsion U, which we assume to be approximately n- 
independent and to interpolate smoothly between the 
limits a(U = 0) = and a(U = oo) = 1. As discussed 
earlier, the magnitude of U depends on the ionic charge 
Z. The relative probabilities of singly and doubly occu- 
pied sites are then 



Pe{n) 



Ps(n) 



Ps{n) +Pd{n) 



P d {n) = 1 - P s (n) 



(50a) 



(50b) 



and the incoherent part of the spectral function in Eq. 
(46) is given by 

A\u) = (1 - z(n))[P s (n)A' h (w) + P d (n)A h (u)]. (51) 

In Figure (10) we show the qualitative behavior of the 
spectral function Eq. (46) versus band filling for one set 
of parameters. As the band filling starts increasing from 
zero, the quasiparticle peak starts decreasing and a peak 
at positive frequencies forms, corresponding to the ex- 
cited states of the doubly occupied ion upon creation of 
electrons in singly occupied ions. Only at higher fillings, 
and particularly beyond half-filling, where the number 
of doubly occupied atoms becomes appreciable, does the 
peak at negative frequencies start to form, corresponding 
to the excited states of the singly-occupied ion resulting 
from destruction of an electron in the correlated ground 
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state of the two-electron atom. The peak at positive ener- 
gies decreases as the full band is approached, because we 
are not considering the possibility of creation of electrons 
in higher energy states of the doubly occupied atom. The 
quasiparticle peak height decreases monotonically as the 
band filling increases, approaching Sf as the band be- 
comes completely full. In the limit a — ► 1, correspond- 
ing to very strong intra-atomic Coulomb repulsion (i.e. 
very large ionic charge Z) the peak at negative energies 
only starts to form for n > 1. Note also, as discussed 
in the previous section, that the relative magnitude of 
the quasiparticle and the incoherent parts will depend 
on the ionic charge Z, with the incoherent contributions 
becoming more important for small Z. Furthermore, the 
energy scale of the incoherent excitations decreses with 
decreasing Z. 



B. Frequency-dependent conductivity 

On general grounds we expect that the frequency- 
dependent conductivity will be given by a <5-function at 
zero frequency broadened to a Drude form by disorder 
and electron-phonon interactions if the spectral function 
is of the ^-function form Eq. (31a), i.e. for the case of 
electrons in an almost empty band. On the other hand, 
as the band becomes full and the single electron spec- 
tral function develops incoherent contributions we expect 
the frequency-dependent conductivity to develop high 
frequency contributions reflecting incoherent absorption 
processes where the atoms i and/or j end up in excited 
states when an electron absorbs a photon and jumps from 
atom i to atom j. A convenient form for the frequency- 
dependent conductivity that reflects this physics is given 
by the form derived within dynamical mean field theory 



/oo />oo 
de / duj'p(e)A(e,oj')A(e,uj + oj') (52) 
-oo J — OO 

/( W )-/( W+W ') 



Here, f(w) is the Fermi function and A(ek,u/) is the single 
electron spectral function. p(e) is the electronic density 
of states for the quasiparticle band. Eq. (52) is exact in 
the limit of infinite dimensions, and is believed to be a 
reasonable approximation for finite dimensions. 

We evaluate Eq. (52) for the phenomenological spec- 
tral function given by Eqs. (46) and (51). For the quasi- 
particle band density of states we choose for simplicity a 
triangular form 



^ = |( 1 + W2) ; -f <£<0 



(53) 



— (1 

IT D/2 



' *><«<7 



with bandwidth D = z(n)Do, and Do the 'bare' band- 
width for electrons in the almost empty band. The 'effec- 
tive bandwidth' D becomes progressively smaller as the 
band filling increases, corresponding to the progressive 
reduction in the quasiparticle hopping amplitude 



t(n) = toz(n) 



(54) 



with to the hopping for the almost empty band. The 
Fermi energy is related to the band occupation by 



e F = —(-i + y/n) ;n < 1 
D 



(55) 



(1 - V2 -n);n> 1 



in this model. The frequency-dependent conductivity has 
three contributions 



C7i (w) = a{(uj) + af(u)) + al(uj). 
At zero temperature, they are given by 



<j\(uf) = ire 2 z 2 (n)p(e F ) — 



i r 



7T w 2 + r 2 



(56) 



(57a) 



ire 2 z(n) f° , 



[A'(u+uj')p(e F +uj') (57b) 



+A'(ui')p(e F +uj + uj')} 



ttp 2 r° 

al(u)) = — dw'A'(u')A'(uj + u'). (57c) 

W J-w 

In Eq. (57a), we have replaced the lo = S- function 
by a Drude-like lorentzian with relaxation time r = 1/T 
as would be expected in the presence of impurities or 
electron-phonon scattering. 

The physical content of the different contributions to 
£7i(o;) is as follows. a\ describes the intra-band opti- 
cal absorption of the quasiparticles, where the atoms re- 
main in their ground state when the electron hops, and 
in particular gives the zero frequency conductivity. Note 
that it is proportional to z[n) = l/m* as expected in 
the Drude model (one of the powers of z(n) is cancelled 
by the density of states factor). a\ has two contribu- 
tions; the first term has a non-zero contribution from the 
incoherent part of the spectral function at positive fre- 
quencies: it corresponds to processes where an electron 
hops to a singly occupied site, and the resulting dou- 
bly occupied atom ends up in an excited state, while the 
atom where the electron hopped from (initially doubly or 
singly-occupied) remains in the ground state; similarly, 
the second term in a\ has a non-zero contribution from 
the incoherent part of the spectral function at negative 
frequencies, corresponding to processes where an elec- 
tron hops from a doubly occupied site and the resulting 
singly occupied atom ends up in an excited state, while 
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the atom the electron hops to (either singly or unoccu- 
pied initially) ends up in the ground state. Finally, af 
describes processes where an electron hops from a doubly 
occupied site to a singly occupied site and the atoms at 
both sites end up in excited final states. 

Figure 11 shows the obtained frequency-dependent 
conductivity for one set of parameters using the phe- 
nomenological spectral function discussed in the previ- 
ous subsection. For low band filling, with the spectral 
function dominated by the quasiparticle part, the opti- 
cal conductivity is largely given by the intra-band Drudc 
part, and the zero frequency conductivity is high. In 
other words, the system is a 'good metal', with small 
resistivity. As the band filling increases, the intra-band 
conductivity first increases as the number of electrons in- 
creases, however at the same time some of the increased 
absorption occurs at higher frequency due to the incoher- 
ent part of the spectral function. At even higher band 
filling where the intra-band conductivity becomes hole- 
like, the Drude part becomes much smaller as does the 
zero-frequency conductivity; at the same time, a large in- 
coherent contribution to the conductivity develops; the 
spectral weight of the incoherent part arises from the re- 
duced intra-band contribution due to the reduced quasi- 
particle weight. The system is now a 'bad metal', with 
high resistivity and largely incoherent optical absorption. 

It should be pointed out that our calculated conduc- 
tivity certainly does not satisfy the 'global' conductivity 
sum rule Jl8[ 

f°° , , , ire 2 n , x 

/ dua^io) = — — (58) 
Jo 2m 

which states that the total conductivity spectral weight 
should increase with the number of electrons in the sys- 
tem (m is the bare electron mass). This is because our 
model does not take into account interband transitions, 
and because it does not allow for absorption processes 
where the atomic final state contains three electrons. 
Nevertheless, the model correctly describes the transfer 
of optical spectral weight from intraband to higher fre- 
quency incoherent absorption as the bandfilling progres- 
sively increases. 

C. Dc conduction 

As a consequence of this physics, it is clear that the 
process of electrical conduction in a hole metal is intrin- 
sically more complicated than that in an electron metal. 
When a hole hops from a site to a neighboring site, there 
is a finite probability amplitude for the atoms to end 
up in excited states rather than the ground state. That 
is what gives rise to the reduced dc conductivity as ex- 
pressed by the enhanced effective mass in the Drudc for- 
mula. Thus the conduction cannot be described as an en- 
tirely coherent process, as in the case of electrons, where 
the phase of the wavefunction is well defined over a mean 



free path. In the case of the hole, only the part of the 
conduction associated with the quasiparticle part of the 
spectral function can be described as coherent. When the 
quasiparticle weight is small, the transport of current is 
largely incoherent. 

It is of course true that in a perfect translationally in- 
variant crystal the many-body eigenstates are also eigen- 
states of the crystal momentum operator po|| , and one 
might argue that such a state with finite k could carry 
current without dissipation. However, any imperfections 
in the crystal would make such arguments invalid, un- 
less the system goes into the superconducting state and 
acquires rigidity. 

One may expect new physical phenomena associated 
with this incoherent transport. If an atom ends up in 
an excited state after the hole hops, it amounts to cre- 
ating an electronic excitation in the system, just as the 
electron-phonon interaction can give rise to creation of a 
phonon . In the electron-phonon case this phonon emis- 
sion gives rise to Joule heating that dissipates part of 
the energy supplied to the system by the external poten- 
tial. Similarly, here the electronic excitation should decay 
through photon emission, with the photon frequency de- 
termined by the electronic excitation energies involved. 
Because of conservation of energy it is clear that such 
photon emission cannot occur in every single hopping 
process for moderate voltages applied. However, one 
would expect emission of non-thermal high frequency ra- 
diation at locations in the solid where translational sym- 
metry is broken, i.e. impurities and boundaries. In par- 
ticular, electroluminescence p2| (at visible and infrared 
frequencies) would be expected to occur in this picture at 
the s&urple-positive electrode boundary, where holes are 
injected into the hole metal sample. The effect should be 
strongest in samples where the holes are highly dressed 
in the normal state, i.e. where the band is close to full 
and the ionic charge Z is small. Such systems should 
also give rise to high T c superconductivity, as discussed 
in the next section. When the system goes superconduct- 
ing transport becomes coherent and the photon emission 
should dissappear. 

Furthermore, the physics of hole conduction discussed 
here suggests that the contact resistance between an 
electron-metal electrode and a highly dressed hole metal 
sample should be non-symmetric. When an electron is 
injected into the hole metal it needs to excite the other 
electron in the orbital, while when an electron is removed 
from the hole metal the remaining electron relaxes to a 
lower energy state. Hence electron injection into the hole 
metal should offer more resistance than electron removal, 
giving rise to a rectifying contact where the resistivity is 
lower at the positive electrode- hole metal contact. 
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VI. HOLE SUPERCONDUCTIVITY 

The physics discussed earlier leads to a mechanism for 
superconductivity [p|-pd| . Consider the low energy effec- 
tive Hamiltonian for quasiparticles Eq. (44): the hopping 
amplitude for a quasiparticle depends on the occupation 
of the sites involved in the hopping process. For zero, 
one and two other electrons in these sites the hopping 
amplitude for an electron is respectively 

t% = Uj (59a) 
t\ 3 = S x Ui (59b) 

hj = Sikj = hi (59c) 

becoming progressively smaller as the occupation in- 
creases, due to the progressive reduction in quasiparti- 
cle weigth. It is convenient to rewrite the quasiparticle 
kinetic energy in terms of hole rather than electron opera- 
tors through a particle-hole transformation (that doesn't 
change the physics). The quasi-hole Hamiltonian is then 

H i h - ~ + ^iji^i,-" + "i,-<r)] 5 L 5 j> + ' 

(60a) 

Atij = t} 3 - t?. = (i- - Ifa = Ttij . (60b) 

Here, the operators are hole rather than electron oper- 
ators. We have also added a term describing density- 
density Coulomb repulsion between holes E/y, of which 
the largest will be the on-site repulsion U =Uu. Further- 
more, we have omitted the term describing hopping of a 
hole in the presence of two other holes, which should be 
negligible for low hole density due to Coulomb repulsion. 

Superconductivity will occur for this model Hamilto- 
nian in the regime of low hole density (band almost full) 
for a range of parameters, driven by lowering of kinetic 
energy when hole carriers pair P,|lC|]. For simplicity we 
now restrict ourselves to the case of nearest neighbor hop- 
ping and interactions only. The hole quasiparticle Hamil- 
tonian is 

H = - [t + &t{n l ,-<y + n 3 ^ a )]{c\ a c ja + h.c.) 

<ij>,a 

i <ij> 

We have omitted the 'tilde' on the hole quasiparticle op- 
erators. This 'correlated hopping' Hamiltonian has been 
extensively studied in recent years by approximate and 
exact techniques, and the reader is referred to the refer- 
ences for detailed information [pp^-|3^|. The condition 



for superconductivity in the limit of low hole concentra- 
tion is m 

K > y/(l + u)(l + w) - 1 (62) 

with K = 2zAt/D, u = U/D, w = zV/D, with z the 
number of nearest neighbors to a site and D the (renor- 
malized) bandwidth. Within BCS theory the same con- 
dition is found, with l/D replaced by <7(£f)j the density 
of states at the Fermi energy, in the definition of the pa- 
rameters in Eq. (62) ||. The reader is referred to the 
literature for discussion of the particular features of the 
superconducting state (e.g. energy dependent gap func- 
tion) and for results for various observables (e.g. density 
and pressure dependence of T c , thermodynamics, tunnel- 
ing, NMR, disorder effects, etc.) that can be understood 
from this low energy effective Hamiltonian ||[l(J . 

Here we wish to focus the discussion on the fundamen- 
tal aspect of the physics related to the high energy de- 
grees of freedom that are not contained in the low energy 
quasiparticle Hamiltonian |ji|,[yj. The key to supercon- 
ductivity lies in the relation between the relative weights 
of the coherent and incoherent parts of the spectral func- 
tion discussed in Sect. Ill and the site occupation, or 
the band filling. As we have seen, as the electronic site 
occupation increases the spectral function becomes more 
incoherent. Consider the Cooper pair wavefunction for a 
hole pair in the full band: 

* = -^EM c alo> (63) 

where |0 > is the full band, and c] CT is a hole creation 
operator. In the limit where the size of the Cooper pair 
(i.e. the superconducting coherence length) goes to in- 
finity the holes become uncorrelated and ~ 1/V~N 
(N =number of sites in the system). For short coherence 
length instead, ~ 0(1) for \i — j\ small. In particu- 
lar, the amplitude for the holes to be on the same site or 
on neighboring sites becomes of order unity rather than 
being negligible. Thus, the spectral function for a hole 
in such a state will have a larger coherent part and a 
smaller incoherent part than the spectral function of a 
hole that is not bound in a pair. In other words, holes 
'undress', i.e. become more coherent, when they pair 
pT[ . Because electrons (^quasiparticles when the band 
is almost empty) are more coherent than holes ^quasi- 
particles when the band is almost full) we may also say 
that holes become more 'like electrons' when they bind 
in a Cooper pair. Naturally, if pairing and superconduc- 
tivity is associated with such quasiparticle 'undressing', 
it will occur when the quasiparticles in the normal state 
are hole-like (band close to full) and will not occur when 
they are electron- like (band close to empty) jj4|. 

Intimately related to the increased coherence of the 
hole bound in a Cooper pair is the fact that the hole has 
an increased hopping amplitude when it hops to or from 
a site where there is another hole, due to its increased 
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quasiparticle weight. When a hole forms a Cooper pair 
the probability for another hole to be nearby increases, 
and the increased hopping amplitude that occurs on the 
average gives rise to a lowering of the hole kinetic energy 
upon pairing, which provides the binding energy of the 
Cooper pair if it can overcome the increased cost in ordi- 
nary Coulomb energy arising from the decreased distance 
between like charges [p5| . 

Finally, consider the process of optical absorption when 
the band is almost full. An incident photon causes an 
electron to hop from a site to a neighboring site. If the 
electron is part of a hole Cooper pair, it is far more likely 
to hop to or from a site that is unoccupied by another 
electron, than if it is uncorrelated. If the electron hops to 
or from an unoccupied site, it has a larger contribution 
from the coherent part of the spectral function and hence 
gives a larger contribution to the 'intraband' part of a\ 
rather than to the high frequency incoherent part. Hence, 
optical spectral weight will shift from the high frequency 
incoherent part of the spectrum to the low frequency 
Drude part when pairing occurs. Because in the super- 
conducting state the lowest frequency part of the Drude 
absorption (below twice the energy gap A) collapses into 
the zero frequency 5— function that determines the pen- 
etration depth, the 8— function acquires an 'anomalous' 
extra contribution transferred from the high energy inco- 
herent optical absorption spectrum rather than from the 
frequency region below 2A, and an 'apparent violation' 
of the conductivity sum rule results |3f| . 

In summary, the superconducting transition driven by 
the physics of electron- hole asymmetry discussed here has 
associated with it these anomalous spectral weight trans- 
fers in one- and two-particle (optical conductivity) spec- 
tral functions, caused by the dependence of the single 
particle spectral function on site occupation discussed in 
Sect. Ill, together with the fact that the site occupation 
changes when pairing occurs. The spectral weight is 
transfered from an energy range that is unrelated to the 
magnitude of the superconducting energy gap A, rather 
it is related to the much larger energy scale of intra- 
atomic electronic excitations. To study quantitatively 
these spectral weight transfer processes it is useful to 
consider model Hamiltonians as discussed in the next sec- 
tion. 



VII. MODEL HAMILTONIANS 

To go beyond the phenomenological model discussed 
in the previous section it is useful to construct model 
Hamiltonians that contain the physics of interest, that 
are amenable to theoretical studies by various techniques. 
We consider here three classes of model Hamiltonians. 



A. Hamiltonians with auxiliary oscillator degrees of 
freedom 

A variety of generalizations of the commonly used Hol- 
stein model describing coupling of a local boson (usu- 
ally phonon) degree of freedom to the local electronic 
charge density can be constructed to contain the essen- 
tial physics discussed here 0. As the simplest example, 
consider an electronic model with a single orbital per site 
and a local boson degree of freedom, with site Hamilto- 
nian 

H^^ + ^Kqi + iU + aqJn^ (64) 

describing the coupling of a local oscillator qi of frequency 
luq = y/K/M to the atomic double occupancy. Complet- 
ing the squares yields 

p^ 1 ct 
Hl= 2M + 2 K( - qi + k 71 ' 1 ^ 2 + Ue ff n ^ n ^ ( 65a ) 



U eff = U- — ^U-u; g 2 (65b) 

so the equilibrium position of the oscillator is q = for 
the empty and singly occupied site, and q = —a/K for 
the doubly occupied site. The reduction from U to U e f / 
due to the change in the oscillator equilibrium position 
parallels the reduction from U to f7 e // in the atomic 
problem due to the modification of the electronic wave 
function in the doubly occupied atom by correlations. 

The states of the site can be written as direct product 
of the electronic state in occupation number representa- 
tion and the boson state. We denote by \n l > the I— th 
excited state of the oscillator when there are n electrons 
at the site (n=0, 1 or 2) , and |n >= |n > the ground 
state. The site ground states are then |0 > |0 >,| t> 
|1 >, | l> |1 >, | H> |2 >. Because the oscillator is not 
affected by single occupancy, \0 l >= |l' >. The matrix 
elements Si defined in Eq. (13c) and excitation energies 
Eq. (13b) are given by 

S m =< 1<|2 >= -^f (66a) 

El = uj Q l (66b) 

and the spectral function for hole creation in the doubly 
occupied site Ah(ui) has the form given by Eq. (14), with 
quasiparticle weight 

z = Sl = e- 9 \ (67) 

The spectral functions for electron creation in the empty 
site and electron destruction in the singly occupied site 
A e i and A' el are single 6— functions, as in the case of the 
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atom, because the oscillator is not affected by single occu- 
pancy. Finally, the spectral function for hole destruction 
(electron creation) in the singly occupied site A' h (uj) is a 
sum of 5— functions like Eq. (21), with 

Si = {-l) l Si (68a) 

Ef = E\ (68b) 

In an atom, the energy level spacing is proportional to 
Z 2 , with Z the nuclear charge. If we assume similarly 

oj = cZ 2 (69a) 

and we also assume 

9 2 = (69b) 

(c, c' constants) then the reduction from U to U e ff in Eq. 
(65b) is independent of Z, which is what is found in real 
atoms Furthermore, as Z decreases g increases and 
the quasiparticle weight Eq. (67) will decrease, similarly 
as found in the atomic problem in Sect. III. For the solid 
we take then as model Hamiltonian 

ff = -Evi c ^+E^ (70) 

ij i 

and argue that this model Hamiltonian describes the 
same physics as the real solid discussed earlier. In partic- 
ular, the low energy effective Hamiltonian is of the same 
form as Eq. (44), and the same processes of spectral 
weight transfer discussed earlier will occur here. Quasi- 
particle operators are easily defined using a generalized 
Lang-Firsov transformation pl] , p3| . 

This Hamiltonian can be studied by standard theoret- 
ical techni ques . Some initial studies have already been 
performed Jll|j8| . Further studies with non-perturbative 
techniques such as dynamical mean field theory [ pl|j36| 
and numerical methods such as density-matrix renormal- 
ization group |57j should be able to determine quantita- 
tively the parameter regime where superconductivity oc- 
curs in this model, as well as provide substantial insights 
into the physics of electron-hole asymmetry and spectral 
weight transfer under discussion here. 

It should be pointed out however that there is a differ- 
ence between this Hamiltonian and the atomic case pre- 
viously discussed. Namely, both the excitation spectrum 
and matrix elements here are the same for hole destruc- 
tion in the doubly occupied site and electron creation in 
the singly occupied site, unlike the case of the real atom. 
However, we do not believe that this difference will qual- 
itatively change the physics. 

It is also possible to extend the Hamiltonian Eq. (64) 
so that it will also describe dressing of electrons in singly 
occupied sites [Q , by including coupling to the electronic 
charge density as in the ordinary Holstein model. As 
long as the coupling to double occupancy is non-zero, or 



through other terms that break electron-hole symmetry 
such as anharmonicity , the essential physics of inter- 
est here, that as the band filling increases the coupling 
strength and hence the dressing of the quasiparticle in- 
creases, will be contained in the Hamiltonian. Such a 
generalized Holstein model would be useful if one is in- 
terested in describing the smaller dressing effects of the 
first electron in a shell due to coupling to electrons in 
inner shells. 



B. Electronic Hamiltonian 

An electronic Hamiltonian without auxiliary degrees of 
freedom that contains the physics of electron-hole asym- 
metry and undressing under consideration here needs to 
have at least two orbitals per site. A simple such Hamil- 
tonian for the site is 

Hi = Un^riii + U'n'^n'^ + Vn l -n! l + en- 

-t'J2(4Aa + h.c.) (71) 

a 

where the primed and unprimed orbitals refer to electrons 
in the upper (lower) energy orbital, assumed orthogonal 
to each other. (An alternative description would assume 
t' = and non-orthogonal orbitals). We assume all pa- 
rameters in Eq. (71) positive, and the ordering 

U' + 2e < V + e < U (72a) 

U,U',V » e >> t'. (72b) 

These conditions ensure that a single electron resides pri- 
marily in the lower level and two electrons reside primar- 
ily in the upper level. Hence this Hamiltonian mimics the 
physics of the Hartree Fock solution of the atomic prob- 
lem discussed previously. The effective on-site repulsion 
defined by Eq. (7) is U e ff — U' + 2e to lowest order in 
t'. 

It is simple to diagonalize the site Hamiltonian ex- 
actly, however for illustration purposes we restrict our- 
selves here to lowest order perturbation theory in t' . The 
ground states of the singly and doubly occupied sites, 
|er » and | fj,» and their energies are, respectively 

\a »= \a > |0 > +<5i|0 > \a >:E = (73a) 

| U»= 10 > | U> W T>U> + U> ' T> ;£ = U' + 2e 

v2 

(73b) 

where the first ket refers to the lower and the second ket 
to the upper orbital, and 

<5i = - (74a) 
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y/2t' 



(74b) 



V -IP - e 

The excited state of the singly occupied atom is 

|cr( 2 ) »= 5i\a > |0 > -|0 > \a >;E = e (75) 
and the excited states of the doubly occupied atom are 

|U (2) » = |T>U> + U>|T> +<5 3 |U>|0> (76a) 
v2 

-<5 2 |0 > | U>;E= V + e 



I U (3) »= I T> 1 1>-| 4> I T>V2; E = V + e (76b) 
I U (4) »= I U> |0 > -<5 3 | T> 1 l> +| |> | T>V2;i? = U 

(76c) 

with £3 = \plt' j (U— V— e). The spectral functions to cre- 
ate an electron in the empty site and to destroy an elec- 
tron in the singly occupied site are single 8— functions, as 
in the atomic case. The 'hole wavefunction' (p a that max- 
imizes the hole quasiparticle weight is given, to this order 
in t', simply by the wavefunction of the upper state. The 
spectral function for hole creation in the doubly occupied 
atom is then of the general form Eq. (14): 



A h {uj) = z5(lv) + Sf <5(w + e) 



1 



V2 



(77a) 



(77b) 



(77c) 



and the spectral function for destroying a hole in the 
singly occupied atom, i.e. creating an electron, is 

Af h (w) = zS(co) + (S' 2 ) 2 5(u -(V-W- e)) 

+ (S' 3 ) 2 S(uj -(U-U 1 - 2c)) (78a) 



S' 2 = -j= - S lS2 



3 v/2 



(78b) 



(78c) 



(78d) 



Hence the site spectral functions are of the same form 
as in the atomic problem discussed previously. The 
'quasiparticle operators' are then or c' ir7 depending 
on whether the site is singly or doubly occupied. We 
may take as lattice Hamiltonian 



ij 
i 

where we have assumed for simplicity the same hopping 
amplitude for the unprimed and primed orbitals. Pro- 
jection of the Hamiltonian Eq. (79) to the lowest site 
energy states yields a low energy effective Hamiltonian 
of the form Eq. (44), which will give rise to supercon- 
ductivity in certain parameter ranges. Furthermore the 
spectral weight transfer processes discussed previously 
should occur for this Hamiltonian. Because it has a rel- 
atively small number of states per site this Hamiltonian 
should be amenable to study particularly by exact diago- 
nalization methods. Parameters in this Hamiltonian can 
be related fairly directly to realistic parameters obtained 
from quantum chemical calculations of real atoms. 



C. Hamiltonians with auxiliary spin degrees of 
freedom 

The physics of interest here was first studied in the 
context of an electronic Hamiltonian with an auxiliary 
spin 1/2 degree of freedom [p|]3§|] to represent the other 
electron(s) in the atom, with site Hamiltonian 

Hi = [F(n it + mi) - ^oK + A<J i + Un^nii (80) 

with V > loq and A << luq,V, and all parameters pos- 
itive. For such parameters, the auxiliary spin in the 
ground state will point predominantly up if — riii — 
and down otherwise. That situation describes then a hole 
Hamiltonian, where the first hole causes a large disrup- 
tion in the 'background' degree of freedom (the spin) and 
the second a small one. We can however also use this 
Hamiltonian with the operators representing electrons if 
we choose the parameter range V < loq < 2V instead, 
so that the first electron will not change the state of the 
spin and the second one will. The lattice Hamiltonian 



H — J]] tij C i(T C j<7 

ij 



5> 



(81) 



also has as low energy effective Hamiltonian the generic 
form Eq. (44) that gives rise to correlated hopping 
p9|| . Exact diagonalization studies of this Hamiltonian 
on small clusters |3q| show that it gives rise to pairing 
in a wide range of parameters. Furthermore, the opti- 
cal conductivity was calculated in a simple case |33|] and 
the results show the spectral weight transfer process dis- 
cussed earlier, as one would expect. 

The site spectral functions are also simple to calculate. 
There is a difference with the atomic case discussed pre- 
viously in that the single electron spectral function for 
the empty site is not a single 5— function here because 
some modification of the background spin degree of free- 
dom occurs even when the first electron is created at 
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the site. Nevertheless, for the parameter range discussed 
the change in the state of the background is much larger 
for the second electron, hence the quasiparticle weight 
is much smaller for a hole than for an electron and the 
same qualitative physics of electron-hole asymmetry and 
undressing occurs. Here again, this Hamiltonian should 
be amenable to detailed study by exact diagonalization. 

Finally, a site Hamiltonian analogous to the oscilla- 
tor Hamiltonian Eq. (64) can also be constructed with 
auxiliary spin degrees of freedom: 

Hi = [Vn^nn - u ]ai + Aai + Un^n^ (82) 

which will contain the physics of interest here for V > luq. 
This Hamiltonian does not 'dress' the single electron as 
Eq. (80) does. A two-parameter version that resembles 
closely the oscillator Hamiltonian would be 

H = uj al + gcv [l - 2n, iT n i Jcr* + Un^n^ (83) 

which is analogous to the oscillator Hamiltonian Eq. (64) 
which in terms of boson creation and annihilation oper- 
ators reads 

Hi = waa\ai + (U + gu (a\ + a^jn^n^. (84) 

For large g, the Hamiltonians Eq. (83) and (84) give rise 
to small quasiparticle weight for the hole. 

In summary, there is a large variety of model Hamilto- 
nians that can be constructed that contain the essential 
physics of interest here, and we discussed some particu- 
lar simple examples. Study of these Hamiltonians should 
yield interesting insights and clarify the physics further. 
However this physics will only exist for certain parame- 
ters in the models, and not for others. In particular, the 
models need to break electron-hole symmetry. Instead, 
if one takes for example the spin Hamiltonian Eq. (80) 
with V = luq, or the Holstein Hamiltonian with only cou- 
pling to the local charge density, the physics of interest is 
lost because the models in those cases are electron-hole 
symmetric. 

VIII. DISCUSSION 

In this paper we have addressed the microscopic rea- 
sons for why holes are not like electrons in condensed 
matter. We have shown that there are qualitative differ- 
ences between electrons and holes in atoms that have fun- 
damental consequences for the physics of the solid state. 
We have proposed that the essential physical difference 
relates to the relation between the coherent (quasiparti- 
cle part) and incoherent parts of the one-electron spec- 
tral function, and the electron versus hole character of 
the quasiparticle. Namely, we have argued that, reflect- 
ing the physics at the atomic level, as the filling of a 
band increases the quasiparticle becomes progressively 
less coherent while at the same time the character of 
the quasiparticle changes from electron-like to hole-like. 



Furthermore, the effective mass of the quasiparticle is 
inversely proportional to the degree of coherence (quasi- 
particle weight) and hence increases as the band filling 
increases. 

This physics naturally leads to a mechanism of super- 
conductivity, since 1) superconductivity involves pairing 
of quasiparticles, and 2) when hole quasiparticlcs pair, 
the local density of electrons decreases. Hence, as holes 
bind in a Cooper pair they become more coherent, and 
more mobile, and more like electrons, and the supercon- 
ducting state is stabilized by the accompanying lowering 
of kinetic energy. As a result the superconducting state is 
more coherent than the normal state, and spectral weight 
at low energies in one- and two-particle spectral functions 
increases while incoherent spectral weight at high ener- 
gies decreases. In other words, quasiparticles 'undress' 
when they pair [fll| . 

The informed reader may wonder about the relation 
between the parameters in the low energy effective Hamil- 
tonian Eq. (44) and the ones calculated in our earlier 
work on hopping amplitudes in the diatomic molecule 
p9[ . For example, Eq, (44) implies that £2/^1 = ti/ta, 
with ti the hopping amplitude for an electron when there 
are i other electrons at the 2 sites involved in the hopping 
process. In fact, the results found in Ref. |ll| do not gen- 
erally obey this relation. The answer is, the calculation in 
Ref. p9] is certainly expected to be more accurate. How- 
ever, the derivation here reflects the essential physics in 
a clearer way and allows to discuss the fundamental pro- 
cesses of spectral weight transfer. The two calculations 
are of course intimately related, in particular recall that 
in the calculations in ref. the difference in hopping 
amplitudes went to zero as the atomic charge became 
large or if one artificially constrained the atomic orbitals 
to not change with occupation, in agreement with the 
physics discussed here. 

One may argue that the difference between electrons 
and holes discussed here is strictly speaking only a 
qualitative difference for the lowest Is band. For higher 
bands, even a single electron in an empty band will 
be partially dressed from interactions with electrons in 
lower bands (not to mention electron-phonon interac- 
tions). While this is true, the considerations in this paper 
certainly imply that within a given band 

z(n)>z(2-n) (85) 

for n < 1, that is, the hole at band filling (2 — n) is always 
more dressed than its counterpart electron at band filling 
n. Hence, the mechanism of superconductivity whereby 
holes in a band pair to 'undress' and become more like 
the electrons in that band should still be operative. 

There is in fact another qualitative difference between 
electrons and holes in energy bands, i.e. between elec- 
trons at the bottom and electrons at the top of the band. 
Namely, the wavefunction for the quasiparticle at the 
bottom of the band has a large amplitude in the region 
between the atoms, thus providing stability to the lattice 
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(bonding state). Instead, the wavefunction for the quasi- 
particle at the top of the band has vanishing amplitude 
in the region between the atoms, as it has highest energy 
and is orthogonal to all the band states below it (anti- 
bonding state). Correspondingly, the amplitude of the 
wavefunction at the atoms is larger for electrons at the 
top than for electrons at the bottom of the band, and as 
a consequence the dressing due to intra-atomic Coulomb 
repulsion also from electrons in other atomic orbitals is 
larger for electrons at the top of the band. 

The phase space of model Hamiltonians that can be 
written down and studied as possible models for con- 
densed matter phenomena is infinite. What makes some 
models more suitable to describe real materials than oth- 
ers? We propose that a useful criterion to validate one 
Hamiltonian over another is whether it can be derived 
from a more fundamental Hamiltonian that describes 
physics on a larger energy scale. In particular: the low 
energy Hamiltonian Eq. (44) (or Eq. (61)) describes 
an increase of low energy spectral weight when carriers 
pair |35|. By itself, it does not contain the high energy 
degrees of freedom where that spectral weight is coming 
from. However, it can be derived from a more funda- 
mental Hamiltonian (e.g. the models discussed in Sect. 
VII) that does contain the high energy degrees of free- 
dom where the spectral weight comes from. In contrast, 
the attractive Hubbard model for example can also de- 
scribe pairing of carriers, and it describes an associated 
decrease in low energy spectral weight (because the ef- 
fective mass increases upon pairing); however, it is not 
clear that it is possible to derive an attractive Hubbard 
model from a more general Hamiltonian that would de- 
scribe a corresponding increase in high energy spectral 
weight. 

As we have discussed, there are a variety of models 
that one can study that contain the physics of interest 
here. Thus, the issue is not the validity of a given model 
Hamiltonian, rather it is the validity of a very general 
paradigm. Namely, that the conventional description of 
energy bands in solids needs to be augmented, beyond 
the knowledge of the e versus k relation and the quasi- 
particle wavefunction, to the description of the behavior 
of the coherent and incoherent parts of the one-particle 
spectral function as function of band filling, and that 
these quantities have the particular dependence on band 
filling discussed here. 

The concepts discussed here are clearly not restricted 
to one particular class of solids but rather are generally 
applicable to all solids, just as the concept of energy 
bands is. Thus, it is unreasonable to expect that the 
mechanism of superconductivity that results from these 
concepts would apply to only one material or any one 
class of materials, or even that it would apply to most 
but not all superconducting materials. In a sense it is 
like superconductivity itself: the collective state is ro- 
bust because it forms a coherent whole, and cannot be 
destroyed by local impurities or imperfections in a given 
solid. Similarly the theory discussed here is a seamless 



web, that cannot be destroyed by individual examples 
that seemingly contradict it: it either applies to all super- 
conductors or to none, and the apparent counterexamples 
(e.g. superconductors that seem to have only electron 
carriers in the normal state, or superconductors that ex- 
hibit seemingly incontrovertible evidence for an electron- 
phonon mechanism) will eventually be explainable within 
the same framework, if the framework is valid. So what 
is the evidence in favor of this framework? 

Most importantly, there is significant empirical evi- 
dence that the presence of hole carriers favors super- 
conductivity p4| . Moreover, the materials with high- 
est T c (high T c cuprates and MgB^) have holes con- 
ducting through negative ion networks (0 = and B~), 
i.e. small ionic charge Z, which favors superconductivity 
within this framework. Moreover, the high T c cuprates 
show clear evidence for 'undressing', both in the normal 
state as holes are added and at fixed hole concentration 
as the system becomes superconducting po|-[45f , as de- 
scribed (in fact predicted HQ) by the present frame- 
work. Moreover, the dependence of T c on doping in 




follows the behavior 
. Moreover, the gen- 



transition metal intermetallics 
predicted by this framework [[47} 
eral empirical relation observed between superconductiv- 
ity and lattice instabilities follows naturally from the fact 
that antibonding states need to be occupied in order to 
have hole carriers. Other experimental evidence in sup- 
port of hole superconductivity is discussed in the refer- 
ences g[H|35]]4|. 

There are many open questions within the framework 
discussed in this paper. A convenient starting point 
for the discussion are the model Hamiltonians discussed 
in Sect. VII. First, are there other, similar or other- 
wise, model Hamiltonians that should be studied, possi- 
bly more appropriate than those? Next one can move 
from there in two opposite directions. In one direc- 
tion, the study of the properties of these Hamiltonians 
by various theoretical techniques. In particular, dynam- 
ical mean field theory [j36| is likely to be a very fruitful 
approach. One should calculate the most basic physical 
properties such as phase diagram, quasiparticle weight, 
spectral functions and various correlation functions. How 
do the basic superconductivity parameters such as coher- 
ence length, superfluid weight, etc. depend on the Hamil- 
tonian parameters? How do the various spectral weight 
transfer processes depend on the parameters? And in 
particular, which observable properties should be calcu- 
lated that can differentiate this mechanism from others? 
In the opposite direction, one would like to connect the 
model Hamiltonians to real materials. What quantities 
should one evaluate in a first-principles calculation of a 
given solid that would be relevant to fix the parameters in 
the model Hamiltonians? Once that question is clarified 
a vast area of study will open up where any given existing 
or proposed material could be tested to see where its pa- 
rameters lie, in particular concerning its superconducting 
properties. 
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APPENDIX A: HYLLERAAS DRITTE 
NAHERUNG 

For the convenience of the reader we summarize here 
the equations that determine the parameters in the 
Hylleraas wavefunction Eq. (5) 0. The energy (in Ry) 
is given by E — Z 2 X, with 



A = — 



L 2 



ANN! 



(Al) 



M 



+ 50ci + 96c 2 + I28cl + 584cic 2 + 1920^ (A2a) 
N = 4 + 35ci + 48c 2 + 96c 2 + 308cic 2 + 576<| (A2b) 



L = 16 - I + (120 - *-) Cl + (192 - ^)c 2 + (A2c) 

(288 - ^)c\ + (1120 - ^p)c lC2 + (2304 - ^)c\ 

Minimization of A yields the coefficients ci and c 2 in the 
wavefunction, and the orbital exponent k is given by 



2M 



(A3) 



APPENDIX B: CALCULATION OF MATRIX 
ELEMENTS S N IN THE VARIOUS 
APPROXIMATIONS 

For the Hartree wavefunction we have for the overlaps 
Eq. (13c) 



S. 



H 



d 3 r(p z ip n (r) = fog, ip n ) = S n (Z, Z) (Bl) 



These integrals are simply evaluated for the wavefunc- 
tions Eq. (26) and yield 



£(-!)<(/ + 2) 



(-)' 



ai 



(B2) 



with ai given by Eq. (26b). 

For the Eckart wave function Eq. (4) we have simply 

E _ S 1 (Z 1 ,Z)S n (Z 2 ,Z) + S 1 (Z 2 ,Z)S n (Z 1 , Z) 

[2(1+£ 1 (Z 1 ,Z 2 )3)]V2 [ 6) 

For the Hylleraas wavefunction, Eq. (6), we use the 
relation 



d 3 r 1 d 3 r 2 F = 2tt 2 / ds du dtu{s 2 - t 2 )F 
Jo Jo Jo 

(B4) 

valid for any symmetric function F(r\,r 2 ), with s,t, u 
given by Eq. (6). The integrals needed are of the form 



I(p, q, r, a, b) 



ds du dtu(s 2 - t 2 )e- as e- bt s p t q u r 

r, Jo Jo 



(B5) 



and yield 

I I 
I (p,q,r,a,b)= bp+q+r+3 [zZ 

q+r+l-j 



J=0 



(q-JV- 



y l 



1 (p + i-l)\(q + r + l- j)\ 



p+i 



P 



(i-iy. 

__ L ___. • , , g!(p + r + l)! - 

and the overlaps are given by 



(B6) 



$n — 



E 



(2Zk) { 



( ^ (f)3) i/ 2 ^i ! 

^ ' b l (Akn) 1 



b e 



1=0 



i=0 v ' 



[1(2 + 1 -i,i,l,a,b) 



1(1 -i,2 + i,l,a,b,) + (-iy +l (I(2 + a, -b) - 



1(1 -i,2 + a, -&))] + ^ [1(2 + l-i,i,2, a, b) 



1(1 -i,2 + i,2,a,b,)- (-iy +l (I(2 + l-i,i,2, a, -b) - 



1(1 -i,2 + i,2, a, -&))] + J [1(2 + l-i,2 + i,l,a,b)- 
1(1 -i,4 + i,l,a,b,) + (-l) l+l (I(2 + 1 - i, 2 + i, 1, a, -b) 



-I(l-i,4: + i,l,a,-b))] 
with 

1 + Ai + A 2 



(B7) 



6 = 



A 2 — Ai 



Ai 



2Zk 



A 2 = 



2kn 



(B8a) 
(B8b) 
(B8c) 
(B8d) 
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FIG. 1. Orbital exponents for the two-electron atom 
versus ionic charge Z for the Hartree (full line), Eckart 
(dot-dashed lines) and Hylleraas (dashed line) wavefunc- 
tions. The dotted line shows the orbital exponent for the 
one-electron atom (= Z). 



FIG. 2. Parameters in the Hylleraas wavefunction Eq. 
(6) describing angular correlations (ci, full line) and radial 
correlations (c2, dashed line) versus ionic charge 



FIG. 3. Effective U (Eq. (7)) versus ionic charge Z for 
the Hartree (full line), Eckart (dot-dashed line) and Hylleraas 
(dashed line) wavefunctions. The dotted line gives the bare U 
(Eq. (8)) for the two electrons in the single-particle orbital. 



18 



FIG. 4. Qualitative depiction of the two-electron state: 
in the Hartree wavefunction, both electrons occupy the same 
expanded orbital; in the Eckart wavefunction, the electrons 
occupy different orbitals; in the Hylleraas wavefunction, the 
amplitude of the wavefunction depends on the relative angu- 
lar and radial coordinates. When removing an electron from 
any of these two-electron wavefunctions, the state of the re- 
maining electron has to change to become the eigenstate of 
the single-electron atom. 

FIG. 5. Orbital exponent of the 'hole wavefunction' ip a , 
defined as the single particle wavefunction that yields max- 
imal quasiparticle weight for the hole spectral function Eq. 
(14), for the various approximations to the two-electron wave- 
function. The orbital exponent of the single electron wave- 
function is also shown (dotted line). 

FIG. 6. Quasiparticle weight for the hole, z, versus ionic 
charge Z, for the various approximate two-electron wavefunc- 
tions. The difference between these curves and 1, the quasi- 
particle weight for the electron, measures the importance of 
electron-hole asymmetry for given ionic charge. 

FIG. 7. Spectral function for hole creation for various val- 
ues of the ionic charge Z. The full lines indicate the magni- 
tude of the coefficients of the (^-functions within the Hylleraas 
wavefunction, the dotted lines and symbols give the results 
for the Hartree wavefunction. The 'width' of the 8— functios 
is arbitrary. As the ionic charge decreases the quasiparti- 
cle peak ((5-function at uj = 0) decreases and the 'incoherent 
parts' (5-functions at negative frequencies) increase. At the 
same time the energy of the incoherent excitations decreases 
in absolute value. 

FIG. 8. Spectral function for hole creation within the 
Hartree approximation for ionic charge Z = 0.4. Here the 
quasiparticle has smaller weight than the incoherent excita- 
tions. 

FIG. 9. Spectral function for hole destruction (electron 
creation) in the singly occupied atom within the Hylleraas 
approximation. Note that the energies where the incoherent 
contributions appear are much lower for hole destruction than 
for hole creation for the same Z. The quasiparticle peak has 
the same height for hole destruction and hole creation. 
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